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1. INTRODUCTION 


Before it is feasible to compute a structural influence coefficient matrix with the 
method of Reference 1 (Appendix B of the Summary Report), the structural properties of 
the wing must be established. In this appendix, an explanation is presented of the method 
used to locate the elastic axis and the method to determine the El and GJ distributions 
along the elastic axes of the wings of the parametric study of Reference 2. The method 
applies exclusively to wings with a 2-spar (front and rear) construction or a single torque- 

box construction. 
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2. SYMBOLS 


The units used for the physical quantities defined in this paper are given both in 
the International System of Units (SI) and the U.S. Customary Units. 


Symbols 

Definition 

Dimension 

AR 

Aspect ratio 


b 

Wing span 

ft (m) 

c 

Average chord 

in (m) 

c r 

Root chord 

ft (m) 

c 

Structural root chord 

in (m) 

s 

E 

Modulus of elasticity 

lb/in 2 (N/m 2 ) 

El 

Bending stiffness 

lb-in 2 (N-m 2 ) 

M A 

J 

G 

Shear modulus / Flexural stress 

lb/in (N/m; 

GJ 

Torsional stiffness 

lb-in 2 (N-m 2 ) 
. 4, 4v 

in (m ) 

I 

Section moment of inertia 

k 

Constant thickness ratio 


l 

Length 

in (m) 

M 

Bending moment 

in-lb (m-N) 

S 

Wing reference area 

ftV) 

t/c 

Thickness ratio 


w 

Load per unit length of span 

lb/in (N/m) 

W 

Total airplane weight 

lb (N) 

\ 

Taper ratio 


A c/4 

Quarter-chord sweep angle 

deg 

a le 

Leading-edge sweep angle 

deg 
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3 . METHOD USED TO LOCATE THE ELASTIC AXIS ON 
TWO-SPAR PLAN FORMS . 


Before the elastic axis can be located, it is necessary to determine where in the 

planform the front and rear spars are located. The result ,s. ° ’ P' anf °™ “Starts 

c- i P K/rsIr-nl finhter winas 1- s reasonable to assume that the elasr c axis srarre 

5fe£ S£Sr. 1 

dependS Fi" ' «7KThr 8 e roun°d"™l« for locating points A and B of the elastic axis 
^ <iT: ^hich the - 

pfanforms with pivots located at 20%, 30%, and 40% half-span (Wings 2 3 4 , ,10, 

the elas^Ec^axl^was as^sumed^toTe ^traTghra^focared^m^way b^twe^nTjeTpars ^ e . 

" ^ 1^1 S^st to^aTa'ighl I'T from ^ centerLe 

f 0 tile pWot/and again a straight line from the pivot to the inboard 50% chord locafon at 

the outboard end of the stroke, see Figure 3. ox , , ■ 

The two planforms with cranked tips (Wings 6 and 8 of Reference :) hoc the r 
.„ tir axis located per Fiqure 2 until the cranked position was reached. At the cranked 
±ll!on Ehe eTaslic a P xU wi assumed to be positioned mid-way between the spars, and vary 
linearly outboard to the wing tip, see Figure 4. 


4. METHOD USED TO ESTABLISH El AND GJ VALUES FOR_ 
TWO-SPAR PLANFORMS 


A typical bending and torsional stiffness (El and GJ respectively) distribution for 
fighter-type wings is presented in Figure 5. It should be emphasized thal t th< e ^' s r b ^ tl ° ns 
of Fiqure 5 apply only to fighter aircraft with a weight of approximately 40,000 lbs and a 
design imit load factor of 7.33. Bending and torsional stiffness distributions similar to that 
of F?aure 5 were generated for all planforms used in the investigation of Reference 1 . 

In using Figure 5 to generate El and GJ plots for different planforms, the stress at 
any point in the new planform must be equal to the stress at corresponding points in the 
planform represented by Figure .5. Development of an expression relating the stress b 
tween two planforms requires the following assumptions: 

1) Stress at any point on the test planform equals the stress at a 
corresponding point on the planform of Figure o. 

2) Thickness ratio (t/c) equals a constant k. 

3) The two wings are made of the same materials. Thus the modulus of 

elasticity (E) and the shear modulus (G) are the same in both wings. 
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Planform Sketch Including Spar Locations 


Figure 1 





Figure 2. Ground Rules for Elastic Axis Location 
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WING NO. 3 OF REF. 1 
SWEPT FORWARD 

PIVOT AT 30% FULLY SWEPT SEMISPAN 

'72° 

70° 

6.05 

0.1 

52.1 84FT (15-906M1 
25.8FT (7.864M) 

‘450 SOFT (41 .8M 2 ) . 

W = 40,000LB (1 .78x1 0 5 N) 
X p 1VOT = ”27.950FT (8.519M) 

W- 4.257FT0.WM) 


SCALE: 1 CM = 2 FT 

(1 CM = .6096 Ml 



2.58FT 

(.7864M) 


26.092FT 

(7.953M) 


Figure 3 


Elastic Axis Location for a Variable Sweep Wing 




WING NO. 6 OF REF I 

TIP CRANKED AT 80% SEMISPAN 

a LE 

t= 45 6 

A c/4 

= 39° 

AR 

= 3.0 

X 

= 0_.25 

c r 

= 20.65FT (6.294M) 

b 

= 38.75FT (11.811 MV 

S w 

= 500 SOFT (46.45m 2 ) 

W 

= 40,0001b (1 .78x10 N) 


5.16FT 
(1 .573M) 





With these assumptions, the stress relationship is derived using the flexural stress equahon: 

„ My (1 


where <r is the flexural stress, M the bending moment y the distance between the neutral 
axis and the outer most fiber of a section and I the section moment of inertia. 

For a cantilever beam with a uniform load, 

...o2 (2) 


where w is the load per unit length of span and c the average chord. The expression for 
the moment becomes: 

m- . (- 

8 $ 

Letting "y", from the flexure formula, equal 1/2 "t" from the thickness ratio, the ex- 
pression for the thickness ratio becomes: 

2y/ c = k ^ 

Letting c = c , the root chord, and dividing both sides of the equation by 2, it follows 


Substituting the expressions for "M" and "y“ into the flexure formula, the flexural stress 

is given by: 9 _ 

Wb ck c 


Equating the stress expressions, it is found that: 


From this relation: 


e, b = 


W A b A : A kc s 


,6S a'a 


<w/s) B b B : B c SE 

(W/ s >A b A S A c s 




W B b B 3 B kc s. 


16S b‘b 


e1 a • 


Equation (?) was used to generate the bending stiffness for the planforms of the parametric 
study re ^ j r J e ^ ( Complexity of the stress formulas for unsymmetric sections, a linear re- 
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El AND GJ 





DISTANCE OUTBOARD OF WING-PIVOT (B.L. 60) ALONG ELASTIC AXIS 
FIGURE 5. Typical Fighter El and GJ Distributions 
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lationship was assumed between the El and GJ values of planforms A and B. 
a ratio was established as follows: 

Let El d be the bending stiffness calculated by Equation (9). Let EI A 
and GJa be bending and torsional stiffness values corresponding to 
Figure 5. For each point on the span represented by Figure 5, new 
GJb values can be calculated for the study planform by using the 
following ratio: 

g j b =<^)gj a 


Thus 


00 ) 


This result indicates that the torsional stiffness is increased by the same percentage that 

the bending stiffness is increased. , _ T _ . ,. , .. e r .. 

The procedure outlined before was used to find the El and GJ distributions for all 

the planforms used in the investigations of Reference 2 . A typical plot for a wing wit 

out pivots is given in Figure 6. 


5. METHOD OF DETERMINING INPUT DATA FOR THE STRUCTURAL INFLUENCE 
" COEFFICIENT SUBROUTINE 


The data used as inputs to the structural program consist of the following: 

1) Coordinates of the end points of the elastic axis segments. 

2) Bending and torsional stiffness, (El and GJ respectively), for each segment 
of the elastic axis. 

3) Assignment of aerodynamic panels to specific end points of the elastic 
axis segments. 

First of all, the elastic axis is broken at the centerline of each chordwise column .of panels, 
the pivot point and the outboard end of the stroke (point A), as shewn in Figure /. Know 
ing the end points of each segment, their respective coordinates can be found, since the 
elastic axis location is known and the y-coordinate of the centerline of each chordw.se 
column of panels can be obtained from the aerodynamic program described in Reference 

(Appendix C of the Summary Report). , . .. . 

Once the coordinates of the elastic axis segments are known, the bending and 
torsional stiffness for each segment can be found with the following procedure. First lo- 
cate the end points of each segment on the corresponding stiffness plot. Figure 8. I hen 
select numerous points between and including the endpoints of each segment For each 
selected point read El and GJ values and take their respective averages. These ElandGJ 
averages represent the El and GJ values for each segment. Table I shows the result of this 

method applied to the planform of Figure 7. . , .. 

The last data needed as inputs to the structural influence coefficient subroutine are 
the elements of the vector matrix IASIGN. The matrix IASI GN Js used to assign aerody- 
namic panels to specific end points of the elastic axis (see also Reference ). In S^ era1 ' 
a panel is assigned to a segment of the elastic axis in a direction perpendicular to the 
elastic axis and is assumed attached to the outboard endpoint of the assigned segment. See 
Fiqures 9 and 10. The panel located below the pivot is assigned to the pivot. The results 
of application of this procedure to the wing in Figure 7 are listed in Table l. 







WING I 
SWEPT 
PIVOT 
a LE = 
A c/4 = 
AR = 
X = 
b = 
c r. = 

Sw 

W 

X PIVOT 

Y PIVOT 

NOTE:- 


NO. 3 OF REF 1 
FORWARD 

AT 30% FULLY SWEPT SEMISPAN 

= 72° 

= 70° 

= 6 . 05 

= 0.1 

= 52.1 84FT (15.906M) 

= 25.8FT (7.864M) 

= 450 SOFT (41 .8m 2 ) 

= 40, 000LB (1 .78x1 0 5 N) 
r = 27.950FT (8.519M) 

i = 4.257FT (1 .297M) 

—m— ELASTIC AXIS SEGMENT 
ENDPOINT 


SCALE: 

1CM = 2FT 


(1 CM = .6096M) 


FRONT SPAR 

! — 


ELASTIC AXIS 



•26.092FT ■ 
• (7.953M) 


Figure 7 Example of Elastic Axis Segments 
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Table 1 . El and GJ Values for the Example Wing of Figure 7 


o -9 
Ib-m x 1 0 


2 _ 6 
N- m x 1 0 


Ib-in^ x 1 0 


1-2 

2- 3 

3- 4 

4- 5 

5- 6 

6- 7 

7- 8 

8- 9 

9- 10 

10 - 11 
11-12 
12-13 


13.19 
13.19 
13.19 
11 .80 
9.10 
6.50 

4.45 
2.55 

1 .45 
0.66 
0.13 
0.060 


37.9 

37.9 

37.9 

33.9 
26.1 
18.65 
12.78 

7.32 
4.17 
1 .894 
0.373 
0.172 


9.59 
7.63 
2.75 
6.10 

7.59 
6.05 
4.00 
2.29 
1 .20 
0.55 
0.085 
0.047 


Table 2. Flemenfr of Matrix "IASIGN 11 for the Example Wing of Figure 7. 


PANELS 


1-10 

11-19 

20 

21-27 

28-30 

31-35 

36-45 

46-55 

56-65 

66-75 

76-85 

86-100 


ASSIGNED TO 
ELASTIC AXIS 
POINT NO. 
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